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In this paper we study regular cosmic string solutions of the non-Abelian Higgs model coupled 
with gravity. In order to develop this analysis, we constructed a set of coupled non-linear differential 
equations. Because there is no closed solution for this set of equations, we solve it numerically. The 
solutions that we are interested in asymptote to a flat space-time with a planar angle deficit. This 
model under consideration present two bosonic sectors, besides the non-Abelian gauge field. The 
two bosonic sectors may present a direct coupling. So, we investigate the relevance of this coupling 
on the system, specifically in the linear energy density of the string and on the planar angle deficit. 
We also analyze the behaviors of these quantities as function of the energy scale where the gauge 
symmetry is spontaneously broken. 
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I. INTRODUCTION 

According to the Big Bang theory, the universe has been expanding and cooling. During its early evolution, 
the universe has been underwent a series of phase transitions characterized by spontaneous symmetry breaking jl|. 
These phase transitions are important, because they provide a mechanism for the formation of topological defects 
such as domain walls, monopoles, cosmic string, among others [2j]. In the eighties of last century the interest to 
cosmic strings has been increase due to the fact that they were considered candidates to provide a mechanism for 
the large-scale structure formation in the universe. Although the recent observations of the Cosmic Microwave 
Background (CMB) radiation have ruled out cosmic strings as seed for the primordial density perturbations, they 
are still candidates to explain a small non-Gaussianity in the cosmic microwave background and the influence on 
the temperature anisotropies. Such effects had origins in the gravitational fields of cosmic strings in motion !3j]. In 
addition, astrophysical and cosmological consequences of formation of strings are emission of gravitational waves and 
high energy cosmic rays such as neutrinos and gamma-rays. These observational data can help to constraint the 
product of the Newton’s constant G to the linear mass density of cosmic string, fj, p}. These effects are generated by 
cosmic strings formed under inflation context which makes the physics of cosmic strings a vast area of interest. 

The first theoretical model associated with strings-like solutions was given by Nielsen and Olesen Q by using an 
Abelian Higgs model lagrangian which presents a spontaneously U (1) gauge symmetry broken. These solutions were 
also named by vortex. The influence of this system on the geometry of the spacetime has been analyzed in [3, H] 
many years ago, by coupling the energy-momentum tensor associated with the Nielsen and Olesen string with the 
Einstein equations. In these papers, static and cylindrically solutions for the matter and gauge fields and also for the 
metric tensor, were investigated numerically by using coupled ordinary non-linear differential equations. A few years 
latter the Abelian Higgs model has also investigated in f| and El- The so-called (p,q) gravitating cosmic string 
have been analyzed numerically in (Tlj. In the latter, static and cylindrically solutions of the system containing two 
different Abelian gauge fields coupled with two bosonic sectors were investigated. 

In this paper we are interested to analyze the cosmic string solution again, however, considering at this time the 
SU( 2) non-Abelian Higgs model coupled to gravity. As it was also shown in to create a non-Abelian topological 
stable string, it is necessary the presence of two bosonic iso-vectors coupled with the gauge field. The interaction 
potential should present terms with fourth powers in these bosonic fields in order to have spontaneously gauge 
symmetry broken. Moreover, it may contain a direct interaction between these bosonic sectors. 

This paper is organized as follows: In section [TT] we present the non-Abelian Higgs model coupled to gravity and 
analyze the conditions that the physical parameters contained in the potential should satisfy in order the system 
provides stable topological solutions. In section Hi Al we present the ansatz for the matter and gauge fields, and for the 
metric tensor. The equations of motion and boundary conditions obeyed by the fields and metric tensor are presented 
in section m In section m we exhibit our numerical results for the behaviors of the fields considering different values 
of the relevant parameters. In this section we present a comparison of this non-Abelian system with the Abelian 
one. Also we investigate the behaviors of the energy density of the string and planar angle deficit, as functions of 
the energy scale where the SU( 2) gauge symmetry is spontaneously broken and the interaction coupling between the 
bosonic sectors. In section [V] we give our conclusions and in the Appendix we present a special solution of the system 
that, although the corresponding configuration is not in according to stable topological configuration, it behaves as a 
BPS solution of the fields equations . 


II. THE MODEL 

The model that we want to study is described by the following action, S: 

S = J (jJ^G R +jCm ') ’ d) 

where R is the Ricci scalar, G denotes the Newton’s constant and C m is the matter Lagrangian density of the 
non-Abelian Higgs model given by 

C m = -±F^F^ a + ±(D^ a ) 2 + ±(D fiX a ) 2 -V(v a ,x a ), a = 1,2,3. (2) 


1 Observational consequences generated by cosmic strings can be found in 11. 
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The field strength tensor is 




[MS 


( 3 ) 


with the SU( 2) gauge potential A*, and e being the gauge coupling constant. The covariant derivatives of the Higgs 
fields are given by [ 6 |], 

D^ a = d^ a + ee abc A b ^ c , (4) 


D^x a = d^x a + ee abc Al x c , ( 5 ) 

where the latin indices denote the internal gauge groups ( a,b = 1,2,3). The interaction potential, y a ), which 

we shall consider is expressed by 

V(<p‘,xl = L [(*>“) 2 -»?] 2 + x [(x “) 2 -oil 2 

+y «*>“) 2 - »?] [( V ) 2 - il ], ( 6 ) 

where the Ai and A 2 are the Higgs fields self-coupling positive constants and A 3 is the coupling constant between both 
bosonic sectors. 771 and 772 parameters correspond to energy scale where the gauge symmetry is broken. The potential 
above has different properties according to the sign of A = A 1 A 2 — A§: 

• For A > 0, the potential has positive value and its minimum is attained for (ip a ) 2 = rfr and (x “) 2 = p 2 - 

• For A < 0, these configuration becomes saddle points and two minima occur for: 


( v >°) 2 = 0 , ( y a ) 2 = r ? 2 + ^ 77 ? 

and 

(x a ) 2 = 0, (v a ) 2 = r,l + ^vl 

M 

The values of the potential for these cases are, respectively, 

4 4 

Vmin = tt-A and V min = —A A. 

4A 2 4Ai 


( 7 ) 

( 8 ) 

( 9 ) 


Both values for V m i n are negatives, since A < 0. 


A. The Ansatz 


The Euler-Lagrange equations for the Higgs and gauge fields, and the Einstein equations for the metric tensor are 
obtained as presented below. 

First let us consider the most general, cylindrically symmetric line element invariant under boosts along z-direction. 
By using cylindrical coordinates, this line element is given by: 


ds 2 = N 2 (p)dt 2 — dp 2 — L 2 (p)d(/) 2 — N 2 (p)dz 2 . 
For this metric, the only non-vanishing components of the Ricci tensor, are: 

NLN" + NN'L' + L(N') 2 

Rtt — Rzz — 7 ? 


( 10 ) 


( 11 ) 


R 


pp — 


2 LN" + NL" 
NL 


( 12 ) 


R<p(/> 


L{2N'L' + NL") 
N 


(13) 
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where the primes denotes derivative with res pec t to p. 

As to the bosonic sectors and gauge fields [12], we use the expressions below: 

/ cos( 0 ) \ 

y? Q 0 ) = f(p) I sin( 0 ) I , 


(14) 


/ - sin( 0 ) \ 

x°(p) =g(p ) I cos( 0 ) j , 


(15) 


and 


A“(p) = 0 


s 

- 

ep J 


3 — 


ep 


9 a , 3 


Wp) = 0, a =1,2, 3. 


We can see that both iso-vector bosonic fields satisfy the orthogonality condition, <p a x a = 0. 


(16) 


(17) 


III. EQUATION OF MOTION 


In order to present the equations of motion in a more appropriate form to apply numerical analyzes, we shall define 
new dimensionless variables and functions as shown below: 


x=y/xiiup, f(p) = r]iX(x), g(p) = r] 1 Y(x) and L(x) = y/ML(p)r] 1 
Therefore, the Lagrangian depends only on dimensionless variables and parameters: 


e rji A 


a = , q = — /3f = v-, * = 1,2,3, 7 = K.rjl and k = 8nG . 

Ai 772 Ai 

It is convenient to use the Einstein equation in the form 

R^v = —k , with T = and p,u = t, x, 0, z 

The energy-momentum tensor associated with the mater field is defined by, 

2 SS 


Tltv yTgSg^ 


9 = det (g^). 


(18) 


(19) 


( 20 ) 


( 21 ) 


Varying the action m with respect to matter fields and metric tensor, we obtain a system of five non-linear coupled 
differential equations. The Euler-Lagrange equations are: 


(. N 2 LX')' 
N 2 L 


= X 


H 2 


x 2 - 1 + ft (v 2 — q 2 ) + -jy 


( 22 ) 


(. N 2 LY')' 
N 2 L 


= Y 


H 2 


{x 2 — 1) + /?2 {y 2 ~ q 2 ) + -jY 


(23) 


) =u{X 2 + Y 2 )H. 


L (N 2 H 
N 2 


As to the Einstein equations (1201) . we obtian: 
(. LNN')' 


= 7 


U -1 x - v - f < y2 - " 2 > 2 - f x - 9 < y2 - « 2 ) 


(24) 


N 2 L 


(25) 
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and 


(. N 2 L')' 
N 2 L 


- H a 
2 aL 2 


(X 2 + F 2 ) |! + I (X 2 - l ) 2 + M (F 2 - q 2 ) 2 + M (x 2 - 1) (F 2 



(26) 


The primes in the above equations denote derivatives with respect to x. 

At this point we would like to notice that the set of differential equation above, Eq.s (|22l) - (l26ll . reduces itself to the 
corresponding one for the Abelian Higgs model 0 , by taking /?2 = Po = 0 and setting one of the Higgs field equal to 
zero. Because one of our objective is to compare our results with the corresponding one for the Abelian case, we shall 
take, when necessary, the bosonic field \ = 0, which, in terms of dimensionless functions, corresponds to take F = 0. 


A. Boundary conditions 


The requirements of regularity at the origin and topologically stable solutions, leads to the following boundary 
conditions for the matter and gauge fields 


H( 0) = 1; H( oo) = 0 


(27) 


X(0) = 0, X(oo) = 1, F (0) = 0, F(oo) = — = q 

m 

and for the metric fields 

N(0) = 1, N\ 0) = 0 , L(0) = 0, L'( 0) = 1 . 
The energy density per unity of lenght is given by: 

£ = J Tt dpdcj) , 


(28) 


(29) 


(30) 


where is the determinant of the (2+1) dimensional metric given in (flOl) by taking dz = 0, and Tq is the 00- 
component of the energy-momentum tensor. Therefore, 


£• = 27T77 i / NL 

Jo 


S + i(^ + n + i(^ + n^ + i (*■-!)■ 


+f + -9 2 ) 2 + f (x»-i)(r»-,>) 


dx . 


(31) 


The and T? components are, respectively, 


T p p = ?7i At 


————- - (X /2 + F' 2 ) + - (X 2 + F 2 ) + - (X 2 - l) 5 
2 aL 2 2 y 2 y ' L 2 4 V ' 


+ M(f 2 -, 2 ) 2 + |(X 2 -1) (F 2 — g 2 ) 


(32) 


T% = vt A 


+ I ( X ,2 + Y' 2 ) - - (X 2 + F 2 ) ^ + - (X 2 - l) 2 
2 aL 2 2 y ’ 2 y ’ L 2 4 V ’ 


+f (F 2 - 9 2 ) 2 + f (X 2 - 1) (F 2 - q 2 ) 


(33) 


and T 2 2 = T*. 
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IV. NUMERICAL RESULTS 

In the following we shall analyze our system. To do that we integrate numerically the equations (1221) - (l26l) with 
boundary conditions m-m by using the ODE solver COLSYS jl5|. Relative errors of the functions are typically 
on the order of 10 -8 to 10 -10 (and sometimes even better). 

As we have said, our objective in this section is to analyze numerically the behavior of the solutions of the non- 
Abelian gravitating string by specifying the set of physical parameters of the system; moreover, we are also interested 
to compare these behaviors with the corresponding one for the Abelian gravitating strings, observing, separately, the 
influence of each system on the geometry of the spacetime. 

Let us first analyze the behaviors of the Higgs and gauge fields besides the metric functions for the non-Abelian 
string system. In order to do that we construct the plots of these fields as functions of the the dimensionless variable, 
x. Our results are exhibited in Fig. |T| In the left panel we present the behaviors of the Higgs, X and Y, and gauge 
fields, H , considering a = 1.0, 7 = 0.6, = 2.0, fa = 1.0 and < 7 = 1 . The behaviors of the metric fields, N and L, 

are exhibited in the right panel. 




FIG. 1: Non-Abelian strings: In the left panel we present the behavior for the Higgs and gauge fields as functions of x considering the 
parameters a = 1.0 7 = 0.6, 02 = 2.0, 03 = 1.0 and q = 1.0. In the right panel we present the behavior of the metric functions. 


Concerning now the behaviors of the matter, gauge fields and the metric functions for the Abelian string, we adopt 
the procedure already explained. So, we plot the fields X and H , besides the metric fields, L and N , as functions of 
x. In the the left panel of Fig. [2] we present the behaviors for the Higgs and gauge fields considering a = 1.0 and 
7 = 0.6. The corresponding behaviors for the metric functions are shown in the right panel. 

The comparison of the behaviors for the metric functions, L and V, as functions of x for both, non-Abelian and 
Abelian strings, are exhibited in Fig. [3l considering a = 1.0, 7 = 0.6, /?2 = 2.0, j3^ = 1.0 and q = 1.0 for the 
non-Abelian string and a = 1.0, 7 = 0.6 for the Abelian string. 

Both, non-Abelian and Abelian string solutions, provide planar angle deficit in the corresponding spacetime. This 
quantity is given by analyzing the slop of L for points very far from the string’s core. The planar angle deficit is given 
by: 


S = 2tt( 1 - L'(oo)) . (34) 

Moreover, the linear energy density is obtained by integrating m- Taking for the parameters the values used in the 
previous plots, we have found for the planar angle deficits and energy densities per unity of ‘inrif, respectively, the 
following results: 


5 na ~ 0.7998 and S A ~ 0.3493 


(35) 
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FIG. 2: Abelian strings: In the left panel we present the behavior for the Higgs and gauge fields as functions of x , considering the parameters 

cx. — 1.0, 7 = 0.6. In the right panel, we present the behavior of the metric functions. 



FIG. 3: This figure shows the comparison between the behavior of the metric functions as functions of x , considering a = 1.0, and 7 = 0.6. For 

the non-Abelian case, we have taken (3 2 = 2 . 0 , ^3 = 1.0 and q — 1 . 0 . 


anc 0 

£na/(^Vi) ~ 1-2856 and ss 1.1646 . (36) 

Another relevant analysis associated with the non-Abelian strings concerns the behavior of their energy densities. 
Two independent analysis have developed by us: The energy density per units of 2nr]l as function of fa and 7 . In the 
left panel of the Fig. [4] we exhibit this behavior as function of fa considering 7 = 0.6. In the right panel we present 
this behavior as function of 7 considering fa = 1.0. For both plots we have taken a = 1.0, fa = 2.0 and q = 1. 

The behaviors of the planar angle deficit in units of 2ir, 5/2it, as function of fa and 7 , are presented in Fig. [5] 
In the left plot we present S/2tt as function of 03 considering 7 = 0.6. In the right plot we present this behavior as 
function of 7 considering /? 3 = 1.0. For both plots we have taken a = 1.0, fa = 2.0 and q = 1. 


The sub-scripts NA and A refer to the non-Abelian and Abelian strings, respectively. 


2 













FIG. 4: 


In the left panel we present the behavior of the energy density per unit of 2nrii as function of 03 considering 7 = 0.6. In the right panel 
we present the behavior of the energy density per unit of 2tvy)^ as function of 7 , considering 03 = 1.0. In both graphs we have taken 

a — 1 . 0 , (3 2 - 2.0 and q — 1. 




FIG. 5: 


In the left plot we present the behavior of the the planar angle deficit in units of 27 t as function of 03 considering 7 = 0.6. In the right 
plot we present this behavior as function of 7 , considering 03 = 1.0. In both plots we have taken a = 1.0, 02 — 2.0 and q = 1. 


The solutions that we are analyzing in this paper, named regular strings, are the ones that present the planar angle 
deficit, 5, smaller than 2n. The planar angle deficit is a measurement of the the modification on the geometry of the 
spacetime caused by gravitational interaction of the system. To obtain regular strings is required an optimal choice 
of the parameters a , 7 , /? 2 , P 3 and q. For fixed values of /? 2 , 03 and q , it is possible to determine regular solutions 
examining the region below the curve in the (a — 7 ) parameter space. However, from equation (1311) . the energy density 
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decreases as a is increased. Moreover, as shown in the Fig. |4(b)| the energy density increases with 7 . Therefore, 
regular strings are obtained until the critical value of 7 is achieved, j cr . This j cr is described by the curve in the 
(a — 7 ) parameter space. If fa = 0.0, j cr is greater than in the case of fa / 0.0 because there is no contribution to 
energy density from the coupled bosonic sectors. This fact is showed in Fig. [Gl where we have considered fa = 2.0 
and <2 = 1 , for fa = 0.0 and 1 . 0 . 



FIG. 6 : This figure shows the the regions in the (a — 7 ) parameter space, which contains regular non-Abelian string solutions. Our plots were 
developed considering two distinct values of ft 3 (see the graph). For both cases, we consider q = 1.0 and ft 2 = 2.0. 


V. CONCLUSION 

In this paper we have investigated gravitating non-Abelian, SU( 2), Higgs model of cosmic strings. In order to have 
non-Abelian stable topological string it is necessary the presence of of two bosonic iso-vectors. As to the interaction 
potential, besides contains forth power terms in the fields, it may present a direct coupling interaction between the 
bosonic sectors. The numerical analysis were developed considering specific values of the physical parameters of the 
model. 

One of our objective was to compare the behaviors of the metric functions corresponding to the non-Abelian and 
Abelian gravitating strings as function of the dimensionless variable, x. These behaviors were exhibited, separately, 
in Fig. |T|and Fig. [2] In Fig [3l we have displayed the behavior of the metric functions, N and L , for both strings. 
We can notice that, for specific values of the parametei0, the non-Abelian string provides a larger planar angle deficit 
than the corresponding Abelian one, although the corresponding energy density is larger. In fact, by our results, the 
planar angle deficit is approximately 2.290 bigger and the corresponding energy density is also bigger, approximately 
1.104. 

As to the non-Abelian string, we have analyzed the behavior of the energy density as function of fa and 7 . In figure 
|4(a)| we can observe that its minimum value occurs for fa = 0, and in the figure [4(b)] we see that the energy density 
increases with 7 . 

The analysis of the planar angle deficit as function of fa and 7 are exhibited in Fig. [5] In the figure [5(a)] we can see 
that the minimum value for this quantity happens for fa = 0 ; of course, in this analysis, it is necessary that 7 > 0 . 

We have also provided in Fig. [Gj a graph that exhibit a region in the (a — 7 ) parameter space where regular strings 
can be formed for fa = 2.0 and q = 1 for two specific values of fa. This region is below the curved line in the plot. 
As we can see, considering fa = 0, for a given a the critical value for 7 which allows regular string is bigger than for 


We have considered a. = 1.0 and 7 = 0.6 for both cases, and ft 2 = 2fts = 2.0 and q = 1 for the non-Abelian string. 


3 
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ft = 1- This fact is consequence of the dependence of the energy density, Eq. ED, with ft. Because the planar angle 
deficit increases with e, the planar angle deficit depends on ft; so vanishing this parameter, to attain 6/2n near unity 
is required larger value of 7 . 

In great part of the analysis developed in this paper we have considered parameters that provide A > 0; however, 
in the Appendix, we analyze a special case of non-Abelian string where this condition is not fulfilled. This solution 
presents some similarity with the BPS solution of the Abelian Higgs string, i.e., it presents the metric function N 
equal to unity, and allows to reduce the set of differential equation to a simpler one. We also analyzed in the Appendix 
the behavior of the BPS Abelian-Higgs string. 
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like to acknowledge the Universidade Federal Rural de Pernambuco. Also the authors would like to acknowledge Betti 
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Appendix A: The Special Case 


As it is well known, the Abelian-Higgs system presents the BPS solution for the fields equation. This solution is 
characterized by the fact of be solution of a set of first order differential equation which minimize the energy density 
of the string. This configuration is achieved by taking a = 2.0 in the equations of motion and assuming N = 1 
0- The BPS limit for the Abelian Higgs model can be obtained from our system of differential equations, by taken 
ft = ft = 0, a = 2.0 and assuming the field Y = 0 and the metric function TV = 1 equal to the unity. In this case 
the set of differential equations reduce to: 


H' = L (A 2 - 1) , L" 


~7 L 


X 2 H 2 

~1X~ 


(* 2 -l ) 2 

2 


(Al) 


and 


X' = X— . (A2) 

1j 

Here, for the non-Abelian string case, we wish to analyze a special solution of the fields equation that present 
characteristic similar to the Abelian BPS one; however, for this case, we have A = A 1 A 2 — A 2 = 0. In order to 
construct this solution we shall substitute a = 2.0, ft = 1.0, TV = 1 and q = 1.0 into the equations (l25l) and (l26l) . So, 
we get: 


H' = L (A 2 + Y 2 ~ 2) 


(A3) 


and 


L" 


~lL 


(A 2 + Y 2 ) 


H 2 

u 


l - (A 2 + Y 2 - 2 ) 2 


(A4) 


Moreover, taking a = 2.0 into Eq. El, two particular independent linear differential equations can be obtained. 
They are: 

X' = X^- and Y' = Y^- . (A5) 

1j 1j 

With these conditions, we can verify that the radial and azimuthal components of the energy-momentum tensor 
vanish: 


T p p = T* = 0 . 


(A 6 ) 


As to the energy density and axial tension, by using equations 

„4 


to 


they can be expressed by: 


T z z = T* = [H( A ' 2 + Y' 2 - 2)]' 


2 L 


(A7) 
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where the derivative is with respect to x. 

Our numerical analysis for both, non-Abelian and Abelian, solutions above mentioned are displayed in Fig. [7] In 
the left panel, we present the behaviors of the Higgs and gauge fields, and the metric functions, for this special case 
for the non-Abelian string, considering a = 2.0, /?2 = P 3 = 1 and q = 1 . In the right panel we exhibit the behaviors for 
the Higgs and gauge fields, and also for the metric functions for the BPS solution of the Abelian strings considering 
a = 2.0. For both cases we took 7 = 0.6. 



FIG. 7: In the left panel we present the behaviors for the Higgs and gauge fields, and for the metric functions considering 7 — 0.6, a — 2.0, 

/3 2 = f3s = 1 and q = 1, for the non-Abelian special case. In the right panel, we present the behaviors of the Higgs and gauge fields, and 
the metric functions for the Abelian BPS string, considering 7 = 0.6 and a. = 2.0. 


Another point that we want to mention is with respect to the energy per unity length. For both configurations, 
the energy density per 27r7^ are equal to unity. As to the planar angle deficit it has been shown in fl 6 | that for the 
Abelian BPS solution, 


6 7 

2 ^ _ 2 ' 

For the special case of non-Abelian string, we can show by combining (1A3I) to (|A5j) that 


(AS) 


L" = -1 [H (X 2 + Y 2 - 2)]' . 


(A9) 


Integrating the above expression we have, 


L'{ 00 ) = 1 — 7 . 


(A10) 


Consequently by (1571) . we obtain, 

— = 7 . (All) 

As we can see, although both solutions presents the same energy density, the planar angle deficit associated with the 
non-Abelian ’’BPS” configuration is twice bigger than the deficit for the Abelian BPS string. The reason for this fact 
is because in the non-Abelian ’’BPS” configuration there are the contribution of two bosonic sectors against just one 
bosonic sector for the Abelian BPS solution in the definition of the planar angle deficit. Finally we want to mention 
that our numerical results have exhibited that: 


• Both energies density per are equal to unity. 

• Both planar angle deficits are in good agreement with the expressions (IA 8 I) and (1A11I) . 
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